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Recently the ATLAS experiment announced a 3 a excess at the Z- peak consisting of 29 
pairs of leptons together with two or more jets, U™ lss > 225 GeV and H t > 600 GeV, 
to be compared with 10.6 ± 3.2 expected lepton pairs in the Standard Model. No excess 
outside the Z-peak was observed. By trying to explain this signal with SUSY we find that 
only relatively light gluinos, nig < 1.2 TeV, together with a heavy neutralino NLSP of 
to* > 400 GeV decaying predominantly to Z-boson plus a light gravitino, such that nearly 
every gluino produces at least one Z-boson in its decay chain, could reproduce the excess. We 
construct an explicit general gauge mediation model able to reproduce the observed signal 
overcoming all the experimental limits. Needless to say, more sophisticated models could 
also reproduce the signal, however, any model would have to exhibit the following features, 
light gluinos, or heavy particles with a strong production cross-section, producing at least 
one Z-boson in its decay chain. The implications of our findings for the Run II at LHC 
with the scaling on the Z peak, as well as for the direct search of gluinos and other SUSY 
particles, are pointed out. 


I. INTRODUCTION 

The discovery by ATLAS [TH3J and CMS [1] experiments at the LHC Collider [5, 6] of a new 
particle with a mass of 125 GeV Ei and with the expected properties of a Higgs boson has 
marked the programme of high-energy physics for the next coming years. On one side, it is 
mandatory to be precise enough in the measurements of the properties of the new scalar particle in 
order to definitively ascertain its nature as the Higgs boson remnant of the electroweak symmetry 
breaking mechanism. On the other, its mass is still compatible with the requirements imposed by 
supersymmetry (SUSY) at the expense of moving the SUSY scale above TeV energies. Combined 
with the current LHC constraints (although model-dependent in most cases) from data analyses 
in the first run of LHC, the search for SUSY effects is becoming more restrictive. Nevertheless, 
the discovery of SUSY would be such an extraordinary event, not only by itself, but for solving 
pending fundamental experimental and theoretical problems in particle physics, that an intense 
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well-motivated experimental programme to search for SUSY effects is of the highest interest. 

Within this scenario, the ATLAS Collaboration has recently presented an intriguing excess 
at the 3 sigma level of e + e~ and n + pairs just at the Z peak [Uj, accompanied by hadronic 
activity and missing transverse energy (MET). With an integrated luminosity of 20.3 fb -1 at a 
p — p CM energy of 8 TeV, the experiment observes a total of 29 pairs of electrons and muons with 
an invariant mass compatible with the Z boson mass, with an expected background of 10.6 ± 3.2 
pairs. No excess over the expected background is observed outside the Z peak 1 . The question 
that immediately arises is whether SUSY, or some other extension of the Standard Model (SM), is 
able to explain that excess of Z + MET events taking into account the current limits on beyond 
the SM physics. A study in those terms within a SUSY framework is presented in this paper. 

As we will show, the observed signal can only be explained if one has a large production cross- 
section of heavy SUSY particles (gluinos or squarks) whose decay chain contains about one Z-boson 
per parent particle. If such an explanation is indeed the answer to the observed excess, our study 
points out the way to confirm it in the Run II of LHC, as well as cosmological implications, in 
particular the particle content of dark matter in the Universe. The resulting scheme of SUSY 
particle mass hierarchy, including charginos and neutralinos, will be apparent. 

II. ATLAS EXCESS IN l+l~ ON THE Z PEAK 

In order to motivate and clarify the assumptions that will be made in the next sections, a 
simplified summary of the mentioned ATLAS search is given here. 

The main focus of the search in [9] are the decays of squarks and gluinos with two leptons 
(electrons or muons), jets and E“ lss in the final state, where the two leptons originate from a 
Z-boson. 

In order to discriminate between SM background events and a possible signal, the following 
requirements are applied: 

• At least two same-flavoured leptons with opposite electric charge are required in each event. 
If more than two leptons are present in the event, the two with the largest values of pt are 
selected. The leading lepton, i.e. the lepton with highest pt , must have pt > 25 GeV, 
whereas the subleading lepton p t can be as low as 10 GeV. Their invariant mass must fall 
within the Z boson mass window, here considered as 81 < m\\ < 101 GeV. 


1 A similar analysis on Z plus E™ lss has been performed by CMS m- However, among other differences, no cut on 
Ht has been applied. No deviation from SM expectations has been observed here. 
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• All events are further required to contain at least two jets with pp > 35 GeV and \r/\ < 2.5, 
to have E!f lss > 225 GeV and Hi > 600 GeV , where Hi is the pi sum over all the jets with 
pi > 35 GeV and \r]\ < 2.5 and the two leading leptons: Hi = Yli^i' 1 +Pp Pt ° n ’ 1 +p!^ pton ’ 2 . 

• Furthermore, the azimuthal angle between each of the two leading jets and £p llss is required 
to be higher than 0.4. 

A great effort has been made to accurately estimate the number of SM events that survive the 
previous selection. An enumeration of these expected SM processes together with some of their 
characteristics follows: 

• The main background, namely it, together with WW, Wt and Z —>• rr, which add up to 
~ 60% of the predicted background, have been estimated using a data-driven method that 
has been thoroughly cross-checked with different techniques. 

• Diboson backgrounds with real Z-boson production (~ 25%) and “rare top” {it + W, it + 
Z, it + WW and t + Z) backgrounds (< 5%) are estimated using MC simulation. These are 
subject to carefully assessed theoretical and experimental uncertainties. 

• Processes with “fake leptons”, i.e. leptons originating from the mis-identification of a jet, 
(~ 10%) are estimated using a data-driven method used regularly in most of ATLAS analyses. 

• Finally, particular care has been taken to suppress the Z/ 7 * + jets background as much as 
possible, given that it could mimic a possible signal (the cut in the azimuthal angle between 
each of the two leading jets and the E™ 1SS has been applied to serve this purpose). Nonethe¬ 
less, a data-driven technique has been used to estimate this small (< 1 %) but important 
background. 

The total number of expected SM model events passing all the requirements is 10.6 ± 3.2 and 
the number of observed data events is 29. This corresponds to a 3.0 a significance. 

III. Z-BOSON PRODUCTION IN THE MSSM 

As explained above, if the observed excess is confirmed, it would clearly point to a new non¬ 
standard process producing additional Z-bosons at LHC energies. Z-bosons are regularly produced 
in the decay chains of most of the SM extensions. Still, this signal would require a significant 
production of Z-bosons without conflicting with all other experimental searches of beyond the SM 
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particles. In fact, using the central value for the expected background and taking into account 
the Z-branching ratio to muon and electron, this would imply that we have produced 273 ± 48 
additional Z-bosons (with 20.3 fb 1 ). Assuming the Z-bosons are produced in the decay-chains of 
beyond the SM particles, Y, produced in the collision, we need to produce at least 273/A f(Y —> Z ) 
Y particles, with A f(Y —> Z ) the average number of Z-bosons produced in the decay of a Y 
particle. On the other hand, as we will see next, the experimental cuts used in the experiment, 
namely 7ij e t s > 2, E™ 1SS > 225 GeV and LLp > 600 GeV , define further the characteristics of the Y 
particle and its decays. 

Now, the following question is whether it is possible to produce these extra Z-bosons with the 
associated event-attributes in some SUSY extension of the SM, while at the same time all the 
constraints imposed by present new-physics searches at LHC and other experiments are satisfied. 
In this section, we will assume that the acceptance of the applied selection, also taking into account 


the reconstruction efficiencies is ideally equal to unity. In Section IV, we will perform a “realistic” 
simulation in a SUSY model using Delphes m to take into account the signal acceptance and 
detector efficiencies. 


A. Production cross-section of MSSM particles 


We need to produce 273 (225 at one-cr) Z-bosons if we want to accommodate the observed excess 
in lepton-pairs. Assuming that R-parity is conserved, supersymmetric particles are produced in 
pairs in processes of the type pp -A YY. Thus, the required cross-section for this process would 
be, 


o(pp -A YY) = 


N ev /J\f(Y -a Z) _ 137(113)/A f{Y -a Z) _ 6.7(5.6) fb 
C “ 20.3 fb” 1 " A f{Y -a Z) 


( 1 ) 


where we take into account that two Y particles are produced in each event. So, if we obtained one 
Z-boson for each Y-particle produced, we would need a production cross-section of (6.7 ± 1.1) fb at 
the LHC with a CM energy of 8 TeV. Then, the first thing we must do is to determine whether it is 
still possible to have these production cross-sections for some supersymmetric particle taking into 
account the constraints from direct searches at LHC. Here, we consider the production cross sections 
of different supersymmetric particles separately to identify the relevant processes. However, in a 


full simulation, as done in Section IV, different sparticles contribute to the final Z plus jets plus 
MET signal and the different contributions must be added. 


Naively, the first option to consider in a hadron collider would be strong production of squarks or 
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gluinos (assuming they produce Z bosons in their decays). However, current experimental searches 
of jets plus missing energy at LHC force the masses of these coloured particles to be high mm- 
Nevertheless, as we will see below, in some cases we can still find cross-sections of the required size. 

Production cross-sections of gluinos and squarks depend only on their masses and are basically 
independent of other MSSM parameters. In the case of gluino and squarks of the first generation, 
the cross-section depends both on the squark and gluino masses due to t-channel contributions, 
but in the case of stops or sbottoms it depends only on the stop or sbottom mass. In Fig. [l] we 
present the production cross-section of gluino pairs and light-flavour squark pairs 2 calculated at 
NLL+NLO with NLL-fast |15H2l| as a function of the gluino or squark mass. Fig. [2] shows the 
squark plus gluino cross-section as a function of the squark or gluino mass with the second mass 
fixed at different values. The different bands in these figures correspond to the cross-section at one- 
a with fixed squark or gluino masses: the blue (dashed) band corresponds to = 1000 GeV, the 
brown (dash-dotted) band to nig t g = 1500 GeV, the orange (dotted) band has nriq j = 2500 GeV and 
the red (solid) band corresponds to decoupled squarks or gluinos. In all these plots, the horizontal 
grey lines represent the required cross-section at one-cr needed to reproduce the referred ATLAS 
excess. 



FIG. 1: On the left, we show the production cross-section of gluino pairs as a function of the gluino mass 
for two fixed values of the first generation squark masses: 1000 GeV (blue/dashed) and decoupled squarks 
(red/solid). On the right, we have the production cross-section of squarks pairs as a function of their mass 
with gluino masses of 1000 GeV (blue/dashed) and 2500 GeV (orange/dotted). 


2 These squark cross-sections are obtained with 5 squark flavours (mass degenerate), i.e. include contributions from 
sbottoms which are treated as a light flavour. 
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FIG. 2: Production cross-section of squark-gluino as a function of the squark mass (left) or the gluino mass 
(right) for fixed values of the associate particle mass: 1000 GeV (blue/dashed) 1500 GeV (brown/dash- 
dotted) and 2500 GeV (orange/dotted). 


As we can see in these figures, the required cross-section is reached only for light gluino and squark 
masses. In the case of gluino production, the needed cross-section is obtained only for mg < 1200 
GeV and favours heavy squark masses. In fact, these gluino masses are in the boundary of the 
allowed region obtained from jets plus missing Et searches at LHC UMl and would contribute 
significantly only if every g produces at least a Z-boson in its decay. For the production of squark 
pairs, present limits are m.q > 1400 GeV for heavy gluinos and mg > 1650 GeV for degenerate 
squarks and gluinos UMl- Under these conditions, a(pp —> qq ) is always well-below the required 
cross-section, even for > 1400 GeV. Here, we do not show the cross-section a(pp —> qq*) as it 
is typically one order of magnitude smaller than a(pp —> qq), and thus irrelevant. 

Another important process is the simultaneous production of squark and gluino shown in Fig. 
[2j However, we see that we would need both squark and gluino to be light, which is not possible if 
we take into account the bounds from LHC searches UMl. In summary, the best option seems 
to be gluino pair production with nig < 1200 GeV with relatively heavy squarks mq > 3000 GeV, 
if we can get at least one Z boson in every gluino decay. 


Yet, we can also consider the strong production of stop pairs, where the current bounds on stop 
masses are much lower, mi > 650 GeV [ 22l [23| . The total tt* production cross-section is shown in 
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FIG. 3: Production cross-section of tt* pairs as a function of the stop mass. The cross-section is basically 
independent of other SUSY masses. 



H t [GeV] 

FIG. 4: Ht distribution in arbitrary units from the decay of a pair of stops of mj. = 743 GeV to jets plus 
U™ ISS after applying the selection Z —> l + l~. 

Fig-d calculated at NLL+NLO with NLL-fast. In this case, we can see that we could reach the 
required cross-section for rrij < 750 GeV which in principle could be achievable in general SUSY 
models (always assuming that every stop produces a Z-boson in its decay). However, the cuts 
Ht > 600 GeV and E™ 1SS > 225 are very restrictive. We can see this in Fig. [4j where we show 
the Ht distribution from the decays of stop pairs with = 750 GeV. As can be seen here, the 
Ht distribution peaks at Ht — 400-500 GeV as a consequence of the relatively small stop mass, 
and only a small fraction of the events are able to overcome the cut on Ht ■ Therefore, we must 
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conclude that it is not possible to generate the required cross-section and fulfil the requirements of 
the observed excess through stop production. 

Apart from the production cross-sections of coloured sparticles, we could still consider the 
weak production of charginos/neutralinos which can be large enough for light gauginos. Tak¬ 
ing into account that the current bounds on chargino and neutralino masses are not very strin¬ 
gent, electroweak production is worth exploring |24| [251 . It is well-known that the largest elec- 
troweak production cross-sections are those corresponding to W and W + W~, in terms of 
gauge eigenstates. This would correspond to X 2 X 1 and X 1 X 1 in terms of mass eigenstates. How¬ 
ever, the observed signal requires Z-boson on-shell, at least two jets, E™ 1SS > 225 GeV and a 
minimum Ht of 600 GeV. As an example, we present, in Figure [5j the production cross-section 
cr(pp —> X 1 X 2 ~t (xi W + ) (xi Z) (xi 33 ) (X? Z 0 )) i n the case where the only restriction imposed 
is a minimum px > 20 GeV for the jets (in blue/dashed), and the same cross-section after applying 
a cut on the hadronic H > 300 GeV . As it can be seen, the cross-section with charginos of 
m + < 350 GeV looks, in principle, able to accommodate the required Z production. However, the 

Xi 

situation changes after we impose the experimental cuts used in the analysis. The required E™ 1SS 
is easily obtained if ?n v 0 > 150 GeV, but the requirement on Ht is very restrictive. In fact, for 
chargino masses m ± ~ 350 GeV, the largest contribution to the production cross-section would 

Vl 

correspond to y/~§ — 700 GeV, which can never produce Ht of 600 GeV and E™ lss > 225 GeV. 
Then, we can expect the cross-section a(jpp —> XiX + ) a t larger values to be strongly reduced. 
This is shown by the red (solid) line in Figure [5j a relatively mild cut on the hadronic Ht reduces 
the cross-section by more than one-order of magnitude. 


Thus we must conclude that, although electroweak production could contribute efficiently to the 
production of additional Z-bosons, these events can not overcome the experimental cuts and can 
not give rise to the observed signal. 

With this channel, we have reviewed all relevant production cross-sections of different super- 
symmetric particles that could potentially explain the signal. The next step would be to calculate 
the average number of Z-bosons per parent particle Y, that we use in Eq. ([I]). 

B. Decay of Supersymmetric particles to Z-bosons 

Z bosons are produced through the decay chains of most MSSM particles, although the number 
of Z-bosons obtained per each supersymmetric particle produced depends on the identity of the 


9 


n 



FIG. 5: Production cross-section of a pair X 2 X 1 as a function of m~± in the CMSSM with mo = 4 TeV, 

Xi 

Ao = 0 and tan/3 = 40 and /i > 0. In blue (dashed), the cross-section with no cut on E™ lss or Ht and in 
red (solid), the cross-section after applying a cut on Hi^ ts > 300 GeV 


supersymmetric particle initially produced and on the supersymmetric spectrum below its mass. 
The main sources of Z-bosons are the decays of neutralinos and charginos and also in some squark 
decays. The couplings of charginos/neutralinos to Z are given by, 

92 


C z xx — 


with 


a Z/J, 

cos % 
yL 


xtr {o$p l + o$p r ) xj + x*v ( of[fp L + ofp R ) y 


1 


O'? = -VnV*, - — Vi 2 V *2 + Sij sin" e w 


32 


o[f = -UtiUj! - -U* 2 U j2 + Sij sin" e w 

= ~l N * N h + 2^A(*4 


o" R = -O 

IJ 1J 


,//L 


( 2 ) 


Therefore, we can obtain Z bosons in the decays of higgsino-like neutralinos and charginos. For 
instance in a usual mSugra spectrum, the second neutralino will only produce Z-bosons through 
its (relatively small) higgsino component while the two heavier neutralinos can be expected to 
produce a sizeable number of Z-bosons. On the other hand, charginos can produce Z-bosons both 
through the wino and from the higgsino component but only in decays of the heavier charginos, as 
the lightest one will only decay to a FF-boson and a neutralino (or lepton-slepton if < m x +). 

In Fig.|6jwe can see the values of Af(x 2 ,3 —> Z) as a function of m x 0 g in an mSugra model. Here, 
we obtain N(x 2 —>• Z) around 0.1, as expected if the higgsino content is relatively small. Then, 
J\f (x3 —> Z) can reach at most 0.45 while the other 50% of the decays go to x ± FF =F . Similarly the 
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heavy charginos produce Z-bosons in their decay as can be seen in Fig. [t} In this case Af(xt %) 
can be 0.3, while we have similar branching ratios to XiW + and Xih- 
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FIG. 6: Average number of Z bosons, J\f{Y —>■ Z), in the decays of X 2 (left) and (right) for a typical 
mSugra spectrum. 
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FIG. 7: Average number of Z bosons, J\f(xt 


Z), in the decays of xt f° r a typical mSugra spectrum 


Besides chargino and neutralino decays, Z-bosons couple also to sfermions through, 


C-Zq = 


-m 


C Tsi - ej sin 2 fa) Z„ + (^ KM )L(^ KM M • 

(3) 


cos Ow 


As we can see, the Z couplings are chirality diagonal and therefore, in decays, they can only enter 
through chirality mixing. Although these couplings could be flavour changing, flavour mixing is 
bounded to be small due to the stringent flavour changing neutral current constraints. Therefore 
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we can expect a sizeable amount of Z-bosons produced only through chirality mixing in third 
generation sfermions in decays like ^ —i► ti + Z, or 62 — > b\ + Z in the large tan (3 regime. 

In addition, we obtain Z-bosons in the decay chains of strongly produced sparticles. We can 
obtain Z-bosons at different steps of the decay chain, either through the couplings of Z to sfermions 
or to charginos/neutralinos that we saw above. For instance, if we produce a pair of gluinos 
a possible decay chain would be g — > t^t hZ t —> Z t xf Z b Z t —y XiW + Z b Z t. 

Therefore taking into account the corresponding branching ratios, this decay chain would contribute 
with 2 x BR(<7 — > x?^ + Z b Z t) to N{g — » Z), the number of Z-bosons produced per g produced. 




1500 2000 2500 3000 


Mj (GeV) 


M- (GeV) 


FIG. 8: Average number of Z bosons, J\f(Y —► Z), in gluino (left) and stop (right) decays for a typical 
mSugra spectrum. 


As we can see in Fig. [8j in a typical MSSM spectrum we obtain at most 0.2 Z-bosons per stop or 
gluino while other squarks produce far fewer Z-bosons per squark. Although Figs. [6}j8] have been 
obtained from a mSugra spectrum, the expected number of Z bosons would be very similar in other 
MSSM versions, as it depends only on the spectrum below the mass of the originally produced Y 
particle and the content of charginos-neutralinos. 

As shown in the previous section, the only possibility to explain the signal is to produce a stop 
or a gluino-pair as the lightest coloured sparticle, being all other squarks much heavier and only 
neutralinos, charginos and possibly some sleptons can be below the gluino or stop mass. Moreover, 
given the size of production cross-sections consistent with the present searches, we need to obtain 
nearly one Z-boson per Y particle produced. Therefore, from the expected numbers of Z-bosons 
that we have seen in this section, we have to conclude that it is not possible to reproduce the 
observed signal in a MSSM with a stable (and light) neutralino. 
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Nevertheless, we can still consider different variations of the MSSM: 

• A first possibility would be to have a light gluino below 1 TeV that can evade the bounds 
from jets plus missing Et if it decays to a sufficiently heavy neutralino LSP in a sort of 
compressed spectrum. Under these conditions the gluino would be abundantly produced at 
LHC and even a small number of Z-bosons per gluino could fulfil the requirements to explain 
the observed signal. However, this would require strongly non-universal gaugino masses and 
very heavy LSP’s and we will not follow this path here. 


• A second option is to consider an MSSM where the lightest neutralino decays to a lighter 
gravitino plus some Z-boson. This is the case in gauge mediated MSSM [26U29] and it could 
be also possible in gravity mediated MSSM if the gravitino is lighter than the neutralino 
which then becomes the NLSP [30] . In this case, the neutralino decays to Z-boson and 
gravitino if it is allowed by phase-space and the branching ratio will depend on the lightest 
neutralino composition. This is the possibility we will explore in the following. 


Thus, we will analyze a situation where the neutralino is the next to lightest supersymmetric 
particle and the LSP is the gravitino. All supersymmetric particles will decay to the lightest 
neutralino which then decays to gravitino plus a photon, a Z-boson or a Higgs. The decay width 
of the lightest neutralino to photon, h or Z plus gravitino [301132] is given by, 

| N\i cos 8\y + IV 12 sin0w| 2 m x 
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1 + 3—£ 
m 


x 


48ttMJ>, 


m 


F(m x ,m d ,m , z ), 


G 


| —N \3 sin o: + Ai 4 cos q| 2 , , 

hG) = -ni+rm- ZZ 2 ~F{m x , rriQ, m h ), 


967rMpj 


rn 


(4) 


G 


with F(x. y, z ) a function of the particle masses non relevant for our discussion that can be obtained 
from m- As we can see, if the lightest neutralino is bino-like, A+i ~ 1, and the mass difference 
between neutralino and gravitino is larger that the Z-mass, the branching ratios are BR(%i —> 
Gj) ~ cos 2 9w — 0.8 and BR(xi G Z) ~ sin 2 9w — 0.2. Similarly if the lightest neutralino 

is wino-like, the branching ratios get exchanged. From this equation we can also see that it 
is possible to get a very large BR to Z bosons as needed to reproduce the observed signal if 
{—Nn sin 9w + A 12 cos 9w) — 1 , but this is only possible if the lightest neutralino has a very large 
wino component. 


Although in gauge mediation models the LSP is always the gravitino, the gaugino masses in 
minimal models are proportional to the gauge couplings and therefore the LSP is mostly bino 
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with small wino and higgsino components. Then, we will have to consider other extensions of the 
gauge mediation idea, like the so-called General Gauge Mediation (GGM) where the gaugino and 
sfermion masses depend on hidden sector current correlators which can be different for different 
gauge groups or particle representations. As we will show in the following section, in these models 
it is possible to obtain general neutralino NLSP as required to explain the excess. 


IV. A POSSIBLE EXPLANATION IN GENERAL GAUGE MEDIATION 


Minimal gauge mediation predicts that all scalar and gaugino masses originate from a single 
scale and powers of the gauge couplings [29] . Recently a model-independent generalization of 
gauge mediation was proposed under the name of General Gauge Mediation E3E3I, where all 
the dependence of soft masses on the hidden sector is encoded in three real and three complex 
parameters obtained from a small set of current-current correlators. In these models the gaugino 
and sfermion masses are given by, 

M r = g 2 M s B l / 2 (0) (5) 

3 

mj = gjY f ( + ^2g 2 C 2 {f\r)M 2 A r , 

r =1 

with 

A = J d y ( 3 C ( i\y) - <>) + Co\y)^j , ( 6 ) 

0\O), Cp r \y) (with p = 0,1/2,1, corresponding to scalar, fermion and vector) are associated 
with the current-current correlators in the hidden sector, £ is a possible Fayet-Illiopoulos term 
(C = 0 in the following), C 2 (f\r) the quadratic Casimirs and M s a characteristic SUSY-breaking 
scale in the hidden sector. 

Having six parameters, B r (0) and A r , to fix the soft masses in the observable sector, it is 
clear now that we have much more freedom in GGM |35H38| and, in particular, we have 


Mi , Mi , M3 

9 / 9 ' 9 ’ 

9 1 g 2 9s 


(7) 


as required to reproduce the observed signal at ATLAS. In particular, we need the NLSP to decay 
to gravitino and a Z-boson with a branching ratio close to one. Fortunately, this is possible in 
GGM as shown in Ref. [38., 39]. 

In this GGM scenario we have used SPheno-3.3.3 @ni an to obtain the full supersymmetric 
spectrum at LHC energies. We define a first parameter point GGM1 with the following parameters, 
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Particle 

9 

Xi 
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X° 

xi 

Xi 

Xi 

G 

Mass (GeV) 

1088.0 

428.4 

431.34 
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1361.0 

4.8 x 10" 9 

Particle 

9L 

9R 

k 
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Ir 

Mass (GeV) 
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2957 
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2952 

2716 

2881 
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Particle 

h 

H 

A 

H+ 





Mass (GeV) 

119.4 

1471 

1471 

1473 






TABLE I: SUSY spectrum in the GGM1 parameter point. 


M s = 400 TeV, b\ /2 = A, = 309 TeV, B 2 1/2 = A 2 = 151 TeV, b\ /2 = 129 TeV, i 3 = 316 TeV 
and tan/3 = 9.8. With these parameters we obtain the spectrum shown in Table [1] With respect 
to this spectrum, some comments are in order: 


1. The two lightest neutralinos and the lightest chargino are very similar in mass, ~ 430 GeV 
and this allows a large neutrino mixing as required. In fact the neutralino mixing matrix is 
given by 


Nij 




/ -0.51 0.85 

0.86 0.51 

-0.015 0.028 
y -0.037 0.065 


-0.076 0.031 ^ 

-0.0024 0.0071 
0.71 0.71 

0.70 -0.71 j 


( 8 ) 


On the other hand the relatively large NLSP mass is needed to overcome the I3™ lss cut. 


2. The gluino is relatively light nig = 1088.0 GeV which allows for a sizeable production cross- 
section and taking into account the squark masses of order ~ 3 TeV. As we can see in Ref. 
[14] . this mass would be allowed at 1 cr in a simplified MSSM with gluino-squark-netralino. 


3. The lightest Higgs mass must reproduce the observed value at LHC of rrih ~ 125 GeV and in 
this spectrum it reaches only 119.4 GeV. This problem (typical in minimal gauge mediation 
models) can be solved either by increasing the stop masses taking a larger A 3 or assuming 
extra operators in the Higgs sector, as the dimension 5 operators proposed by Dine, Seiberg 
and Thomas [42] , Here, we assume that this problem is solved by one of these mechanisms, 
given it does not affect the observed phenomenology on the Z-peak 
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FIG. 9: Ht (left) and E™ lss (right) distributions after applying the selection detailed in section [TTj except 
for the cuts on H t and E™ lss . The GGM1 point is represented by the dashed red line and the GGM2 point, 
by the solid blue line. 



FIG. 10: E™ lss distribution corresponding to the GGM1 point from strong production (solid blue line) and 
electroweak production (dashed red line) with the same selection as in Fig. [9j after applying the Ht > 
600 GeV cut. 

Under these conditions, the lightest neutralino width is T Vl = 4.097 x 10~ 10 GeV and the decay 
branching ratios are BR(x° —> Gy) = 1.14 x 10 3 , BR(%° -ft GZ) = 0.997 and BR(%° —* Gh ) = 
1.35 x 10” 3 . Therefore, gluinos are produced at LHC with a cross-section of 8.4 ± 1.6 fir at 
NLL+NLO and after going through different decay chains all of them produce a Z-boson plus a 
gravitino. In this case, the strong production represents approximately 20% of the total production 
cross-section. 

We simulate the production of supersymmetric particles at LHC at 8 TeV (LHC8) with this 
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spectrum using Pythia 8.1 |43| with Prospino2 [13 H E] K-factors and the response of the 
ATLAS detector using Delphes m- The selection of events for this study is done as close as 
possible to that performed in ATLAS |Jj. The dashed red line in Fig. [9] and Fig. 10 shows the 
Hi and E!f lss distributions respectively for the GGM1 point after applying all the selection detailed 
in section |TTj except for the Hi and E™ 1SS cuts. In the H x distribution, we can distinguish the two 
peaks corresponding to electroweak production at low H x values and gluino production at higher 
Hi . From here, we can expect that the cut on Hi will eliminate most of the electroweak production 


but not the gluino production. This can be seen in Fig. 10, where the I?™ 1SS distribution is presented 
separately for strongly produced events (solid blue line) and for the electroweak component of the 
production (dashed red line) for the same selection as in Fig. [9] but after applying the Hi cut, i.e. 
the final selection except for the cut on E!f lss . The electroweak component is significantly reduced 
by the Hi cut while mainly only events coming from strong production survive the cut, as expected. 
In fact, in this simulation of point GGM1, 99% of the gluino points and only 11% of the electroweak 
points have survived the Hi cut. We see that the peak in the Elf lss distribution is approximately 
at m x o/2, and, for m x o = 425 GeV, a reasonable fraction of the events will survive the E™ 1SS cut 
at 225 GeV. In the simulation, 65% of the gluino point and 53% of the electroweak points survive 
this cut. However, due the relatively small production cross-section, the final number of events 
is small. In this simulation and after applying all relevant experimental cuts, an expected signal 
of 6.34 ± 1.02 lepton pairs is found, to be compared with the observed excess of 19.4 ± 3.2. This 
number of surviving events was obtained at NLO with Pythia and Prospino2 but, unfortunately, 
it is still too low to explain the observations. 

Trying to obtain a model able to account for the excess, we consider a second point in our 
GGM scenario with a lighter gluino. The GGM2 point is obtained with M s = 400 TeV, B ^ 2 = 
A\ = 309 TeV, b)J 2 = A 2 = 150 TeV, b \[ 2 = 110 TeV, A 3 = 270 TeV and tan/3 = 9.8. With 
these parameters we obtain the spectrum shown in Table [0} The neutralino mixing matrix in point 
GGM2 is similar to the corresponding mixing matrix in GGM1, and the BR(%5 —> GZ ) = 0.94. 
However, gluino is now much lighter and the gluino-gluino cross-section is now (41.6 ± 7.5) with 
Prospino at NLL+NLO, thus we can expect many more gluino pairs to be produced and a larger 
contribution in the final selection for this point. 

The simulation for this GGM2 point is presented by the solid blue line in Fig. [9| The Hi 
distribution peaks at slightly lower values than in the case of GGM1, due to the slightly lower 
gluino mass, but it is still enough to overcome the Hi cut at 600 GeV. The E!f lss distributions are 
similar for both GGM points, due to the very similar neutralino masses in both cases. However, 
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Xi 

X° 
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X° 

xt 

xi 

G 

Mass (GeV) 
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1117.1 
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Particle 

9L 

QR 
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^2 



II 
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Mass (GeV) 
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2400 
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2250 
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5360 

Particle 

h 

H 

A 

H+ 





Mass (GeV) 

118.1 

1250 

1250 

1253 






TABLE II: SUSY spectrum for the point GGM2. 


in the case of GGM2 the strong production cross-section is larger and much more important in 
relation with the electroweak production: for the GGM2 point, the strong production represents 


~ 55% of the total cross-section. As seen in Fig. 11, after applying the selection we obtain an 
expected signal of 28.0 ± 4.7 events, slightly over the excess reported by ATLAS, showing that a 
signal point defined along these characteristics can be able to reproduce the excess. 



FIG. 11: E™ lss distribution corresponding to the GGM2 point from strong production (solid blue line) and 
electroweak production (dashed red line) with the same selection as in Fig. [9j after applying the H t > 
600 GeV cut. 


We have to emphasize here that it is not difficult to obtain the observed excess for light gluino 
masses, and a gluino mass in between the two presented examples, rrig E [900,1100] GeV, could 
reproduce the observed signal. However, these points may be in conflict with direct searches of 
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TABLE III: SUSY spectrum for the point GGM3. 


jets plus E™ lss mm- There is a tension between this excess and the bounds from gluino searches 
in jets plus E!f lss . To quantify more accurately this tension, we use the program Checkmate 
that allows to compare the different points of the model with various experimental analyses 
determining whether the point is excluded or not at 95% C.L.. We applied the constraints on jets 
plus E™ iss from ATLAS studies: ATLAS.1308.1841 g7], ATLAS.1308.2631 gS], ATLAS.1407.0583 
[23] and ATLAS.1405.7875 [l3j. In fact, we obtain that this point is indeed excluded by the analysis 
ATLAS. 1405.7875, in the signal region with 6 (or more) jets, with an r = — = 2.0, where S 
is the total number of expected signal events, AS is the total lcr uncertainty on this number and 
Sql is the experimentally measured 95% confidence limit on signal events. 

Thus, we have to find another parameter point able to provide the signal, but satisfying all 
the present constraints on jets plus E™ 1SS . The main difficulty in the GGM scenario is that each 
gluino decay produces typically four or more jets with large pt and E™ ,ss . Therefore, these points 
clash with constraints from observables with a large number of jets (plus E™ 1SS ), which have low 
backgrounds in the Standard Model. On the contrary, constraints with fewer jets are easier to 
satisfy because of larger backgrounds. Then, our strategy will be to compress sufficiently the 
spectrum to reduce the number of (observable) jets. This can be done by increasing the lightest 
neutralino mass closer to the gluino mass, so that the jets in the decays g — > xi + jets have 
smaller px • With this goal, we construct a third GGM point. The GGM3 point is obtained with 
M s = 4160 TeV, B \ /2 = A t = 662.5 TeV, B ^ /2 = A 2 = 344 TeY, b \ /2 = 117 TeY, A 3 = 330 TeV 


and tan j3 = 34.4. With these parameters we obtain the spectrum shown in Table |III[ For this 
point, the BR(x) ) —> GZ) = 0.98 and the gluino production cross section is (22.8 ± 3.3) fb with 
Prospino at NLL+NLO. In this case, the signal is lower than in GGM2 but we still get (13.1 ± 2.2) 
events, well above the experimentally measured 95% confidence limit on signal events. 


Concerning the experimental searches on jets plus MET, according to Checkmate this point is 
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allowed with an r = 1.0 with the best signal region being 2 jets plus E!f lss from ATLAS.1405.7875. 
Moreover, in the equivalent search by CMS of two leptons, jets, and E ™ 1SS , CMS-SUS-14-014 [TO] , 
we obtain r = 1.2 in Checkmate. This result shows some tension with the ATLAS possitive signal 
from the bin with large E™ 1SS in the CMS analysis, but still consistent with observations slightly 
above 95 % C.L.. 



FIG. 12: E™ lss distribution corresponding to the GGM3 point from strong production with the same selection 
as in Fig. [9j before (red-dashed) and after (solid-blue) applying the H T > 600 GeV cut. Notice that 
electroweak production is negligible for the neutralino masses in GGM3. 


Therefore, we have proved that it is indeed possible to construct a Supersymmetric model that 
accommodates the observed excess of lepton pairs on the Z-peak. The simulations presented here 
are only a proof of existence and the final model may be very different. Nevertheless, this model 
will have to share the main features of the examples that we presented here. 

V. PROSPECTS FOR SUSY AT LHC13 

As we have shown in the previous sections, the excess observed in ATLAS, if due to SUSY, 
requires a gluino of a mass ~ 1 TeV producing nearly one Z-boson per gluino in its decay. This 
scenario would also require relatively heavy squarks of the first generation with rriq > 2.5 TeV. If 
this is indeed the correct explanation to the observed excess, such light gluinos would be abundantly 
produced at Run II in LHC together with other SUSY particles. Therefore, the results obtained 
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at pp collisions at 13 TeV will immediately confirm or reject this supersymmetric explanation of 
the ATLAS excess. 

Using NLL-Fast with the spectrum of point GGM1, the gluino pair production cross-section at 
8 TeV was a (pp — > 9ff)LHcl NL ° = 7.6T1.3 fb at NLL+NLO. Similarly the production cross-section 
at 13 TeV would be a(pp — > 55)^013^° = 1504=16 fb, that is, we would expect to produce 20 times 
more gluinos at LHC 13 for the same integrated luminosity. Repeating the same exercise with point 
GGM3, we have a gluino pair production cross-section at 8 TeV of cr(pp —> = 22.8 ± 

3.3 fb while the production cross-section at 13 TeV would be cr(pp —»• SsOlhcis^ 0 = 844 ± 44 fb. 
In this point, the cross section increases by a factor ~ 15 in going from 8 to 13 TeV center of mass 
energy. In any case, for both points this would have unambiguous signatures, both on the .Z-peak 
with a scaling of the signal found at LHC8 and in direct searches for gluinos using jets plus missing 
Et in the extension of the analysis in mm- 

To confirm that it is indeed SUSY behind the excess found at the .Z-peak, we should search 
for other sparticles at LHC13. As shown in Section III A[ squarks of the first generation are 
preferred be heavy to increase the gluino production cross-section, namely above 2.5 TeV. Then 
their production cross-section at LHC 13 would be around 1 fb which would make direct detection 
very challenging. Charginos and neutralinos are expected to be bellow the gluino mass and in some 
cases could be abundantly produced at LHC. Although our analysis does not fix the masses of xf 
and Xi, 2 ) they are expected to be rather heavy m x o > 300 GeV and degenerate. In some cases, 
the production cross-section could be large. For instance cr(pp —>• Xi~Xi)lhcs = 60.2 =t 0.9 fb and 
c t(pp —> Xi"Xi)lhci 3 = 154 ± 5 fb calculated with Prospino for m x o = 313 GeV and r/y = 314 
GeV (with large W component in Xi)- Notice that, as explained above, electroweak production 
is also large for these points at 8 TeV, but is eliminated by the cuts on Hr and number of jets. 
Thus, a large electroweak production could be expected at LHC 13 and dedicated searches should 
be encouraged, specially taking into account the requirement of a large production of Z-bosons 
in chargino and neutralino decays. Finally, the signal does not constrain the masses of sleptons 
or third generation squarks and, a priory, we can not make any prediction on their production at 
LHC13. 

Before closing, we should comment on the nature of dark matter in our scenario. As we have 
seen the signal seems to prefer a non-stable neutralino decaying to a very light gravitino and a Z- 
boson. Under these conditions the neutralino mass has no relation with the dark matter abundance 
of the universe and its role as dark matter compon ent would be played by the gravitino. The 
gravitino mass is not bounded by the observed signal but regardless of its exact mass, unfortunately 
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no signal of dark matter is to be expected in direct search experiments. 

VI. CONCLUSIONS 

Recently the ATLAS experiment announced a 3 a excess at the Z-peak consistent of 29 pairs 
of leptons observed to be compared with 10.6 ± 3.2 expected lepton pairs. No excess outside the 
Z-peak was observed. By trying to explain this signal with SUSY we found that only relatively 
light gluinos, rrig <1.2 TeV, together with a heavy neutralino NLSP of m % > 400 GeV decaying 
predominantly to Z-boson plus a light gravitino, such that nearly every gluino produces at least 
one Z-boson in its decay chain, could do it. 

Unfortunately, this is not possible withing minimal SUSY models, as rnSugra, minimal gauge 
mediation or anomaly mediation. The requirement of a neutralino NLSP decaying to Z plus 
gravitino points to models of General Gauge mediation as the simplest possibility. We have shown 
that a model of this class is able to reproduce the observed signal overcoming all the experimental 
cuts. Needless to say, more sophisticated models could also reproduce the signal, however, they 
will ALWAYS share the above mentioned features, i.e. light gluinos (or heavy particles with a 
strong production cross-section) with an effective A f(g -> Z) ~ 1. 
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